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THEORY*  OF  WING  OF  CIRCULAR  FLAK  FORM* 

By  N.  E.  Kochin 

A theory  is  developed  for  a wing  of  circular  plan  form.  The  dis- 
tribution of  the  bound  vortices  along  the  surface  of  the  wing  is  con- 
sidered in  this  theory,  which  has  already  been  applied  in  a number  of 
papers.  In  particular,  the  case  of  the  circular  wing  has  been  examined 
by  Kinner  in  reference  1.  ~ 

A second  method  is  considered  herein  which  permits'  obtaining  an 
expression  in  closed  form  for  the  general  solution  of  this  problem. 

The  wing  is  assumed  infinitely  thin  and  slightly  cambered  and  the  problem 
is  linearized  in  the  usual  manner. 

Comparison  of  the  results  of  the  proposed  theory  with  the  results  of 
the  usual  theory  of  a wing  of  finite  span,  shows  large  divergences, 
which  indicate  the  inadequacy  of  the  usual  theory  of  the  case  under 
consideration.  For  the  wings  generally  employed  in  practice,  which 
have  a considerably  greater  aspect  ratloj  a more  favorable  relation 
should  be  obtained  between  the  results  of  the  usual  and  the  more  accurate 
theory. 


1.  Statement  of  the  Problem 

The  for’-'ard  rectilinear  motion  of  a circular  wing  with  constant 
velocity  c is  considered.  A right-hand  system  of  rectangular 
coordinates  Oxyz  is  used  and  the  x-axis  is  taken  in  the  direction 
of  motion  of  the  wing.  The  wing  is  assumed  thin  with  a slight  camber 
and  has  as  its  projection  on  the  xy- plane  a circle  ABCD  of  radius  a 
with  center  at  the  origin  of  the  coordinates  (fig.  2,  in  which  a section 
of  the  wing  in  the  xz-plane  is  also  shown). 

Let 

z = !;(x,y)  (1.1) 

represent  the  equation  of  the  surface  of  the  wing,  where  the  ratio  £/a 
as  well  as  the  derivatives  d£/dx  and  d£/dy  are  assumed  to  be  small 
magnitudes. 


*"Teoriya  kryla  konechnogo  razmakha  krugovoi  formy  v plane." 
Prikladnaya  Matematika  i Mekhanika,  Vol.  IV.  No.  1,  1940,  pp.  3-32. 
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The  coordinate  axes  are  assumed  to  be  Immovably  attached  to  the  wing. 
The  fluid  Is  considered  Incompressible  and  the  motion  nonvortical,  steady , 
and  with  no  acting  external  forces.  The  velocity  potential  of  the 
absolute  motion  of  the  fluid  will  be  denoted  by  9(x,y,z)  so  that  the 
projection  of  the  absolute  velocity  of  a particle  of  the  fluid  is  deter- 
mined by  the  formulas 


v „ V » v -2* 

* ax’  r dy'  * bz 


(1.2) 


The  equation  of  continuity 


dv  bv  bv, 


shovb  that  the  function  9 must  satisfy  the  Laplace  equation 

b2*  b2*  b2* 

n ay2  dz2 


(1.3) 


At  the  leading  edge  of  the  wing  the  velocity  of  the  fluid  particles 
is  assumed  to  become  infinite  to  the  order  of  l/V^  yhere  5 is  the 
distance  of  the  particle  to  the  leading  edge;  at  the  trailing  edge  the 
velocity  is  assumed  finite.  From  the  trailing  edge  of  the  wing  a 
surface  of  discontinuity  is  passed  off  on  which  the  function  f suffers 
a discontinuity.  The  function  <r(x,y,z)  and  all  its  derivatives  over 
the  entire  space  bounded  by  the  said  cuiface  of  discontinuity  and  the 
surface  of  the  wing  are  ^continuous. 

The  problem  is  linearized  in  the  following  manner.  The  function  9 


is  assumed  to  suffer  a discontinuity  on  an  infinite  half-strip  J located 

in  the  xy-plane  in  the  direction  of  the  negative  x-axis  from  the  rear 
semicircumference  BCD  of  the  circle  S to  infinity.  In  the  same 
manner,  the  condition  on  the  surface  of  the  wing  is  replaced  by  the 
condition  on  the  surface  of  the  circle  "S  located  in  "the  xy-plane 
and  in  this  way  the  function  <p(x,y,z)  is  assumed  to  be  regular  in  the 
region  obtained  by  cutting  the  inf i nite "half- strip  Z ~and  the  circle  S 
from  the  entire  infinite  space. 
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The  boundary  condition  must  be  satisfied  on  the  surface  of  the  wing. 


~ = c cos(n,x) 
on 


(1.4) 
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where  n is  the  direction  of  the  normal  to  the  surface  of  the  wing. 

The  direction  of  this  normal,  because  of  the  assumption  of  small  curva- 
ture of  the  wing,  differs  little  from  the  direction  of  the  z-axis.  If 
small  terms  of  the  second  order  are  rejected  according  to  the  formula 


cos(n,x 


in  place  of  equation  (1.4), 


(1.5) 


dz 


This  condition  must  be  satisfied  on  the  surface  of  the  wing,  but 
it  is  assumed  satisfied  on  the  surface  of  the  circle  S>  that  is,  for 
z as  Oj  this  again  reduces  to  the  rejection  of  small  teres  of  the  second 
order  by  comparison  with  those  of  the  first  order. 


The  boundary  condition  is  obtained: 


(&\  - - c for  x2  + y2<  a2  (1.6) 

'°Z'Z=0 

which  must  be  satisfied  on  both  the  upper  and  lower  sides  of  the 
circle  S. 

The  boundary  conditions  arc  set  up  which  must  be  satisfied  on  the 

surface  of  discontinuity  Z . On  the  surface  of  discontinuity  at  the 

trailing  edge  of  the  wing,  the  kinematic  condition  expresses  the  con- 
tinuity of  the  normal  component  of  the  velocity,  that  is,  the  magnitude 
bvfbn.  must  remain  continuous  in  passing  through  the  surface  of  dis- 
continuity. Since  on  the  surface  of  discontinuity  the  direction  of  the 
normal  differs  little  from  the  direction  of  the  z-axls,  transfer  of 
the  condition  on  the  surface  of  discontinuity  to  the  half-strip  Z , . 
gives  the  equation 

(1)^+0  (ILo  for  M < a}  + y2  > a2j  x « 0 <x-7> 

which  expresses  the  continuity  of  c^/^z  in  passing  through  the  surface 
of  discontinuity  Z . 
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The  dynamical  condition  expressing  the  continuity  of  the  pressure 
in  passing  through  the  surface  of  discontinuity  £ is  considered. 

In  order  to  determine  the  pressure  p,  the  formula  of  Bernoulli  is 
applied  to  the  steady  flow  about  a wing  obtained  by  superposing  on  the 
flow  considered,  a uniform  flow  with  velocity  c in  the  direction  of 
the  negative  x-axis.  In  this  steady  flow  the  velocity  projections  are 
determined  by  the  equations 


*9 


v*-"c+£}  yy 


b<P 

by 


bj 
2 " bz 


and  therefore  the  formula  of  Bernoulli  reduces  to  the  -fora 

P - - 2 [(-  C \ I)'  + (I)'  + (If)2]  + C°n8taUt  (1*8) 

Rejection  of  Bmall  terms  of  the  second  order  results  in 


bcp 

P^^PC- 


(1.9) 


where  Pq  is  the  value  of  the  pressure  at  infinity. 

Since  the  pressure  must  remain  continuous  in  passing  through  the 
surface  of  discontinuity  at  the  trailing  edge  of  the  wing,  the  equation 
obtained  shows  that  bp/bx  does  not  suffer  a discontinuity  on  the 

surface  of  discontinuity.  Transfer  of  this  condition  to  the  surface  Z 
yields  the  condition 


|y|  < a;  x2  + y2  > a2;  x < 0 (1.10) 


which  expresses  the  continuity  of  bp/bx  in  passing  through  £ . 


The  function  <p  suffers  a discontinuity  on  the  surfaces  S 

and  £ , which  means  that  along  the  surfaces  S and  Z , surface 

vortices  are  located  as  shown  in  figure  2.  The  direction  of  such  a 
surface  vortex  is  perpendicular  to  the  direction  of  the  relative  velocity 
vector  of  two  particles  of  the  fluid  adjacent  to  the  surface  of  dis- 
continuity on  its  two  sides.  In  particular,  on  the  surface  on 

account  of  equation  (1.10),  only  bcp/by  suffers  a discontinuity  and 

therefore  the  vortex  lines  on  Z are  directed  parallel  to  the  x-axis 
as  shown  in  figure  2. 
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Since  all  the  vortices  lie  in  the  ~xy- plane,  at  tvo  points  synmetri- 
cal  vith  respect  to  the  xy- plane,  the  values  of  d?/dz  will  be  the  same, 
whereas  the  values  of  &p/ox  and  d<p/$y  will  differ  only  in  sign. 

It  may  therefore  be  assumed  that 

9(x,y,-z)  - - *(x>y>*)  (1-11) 

Assuming  in  particular  z e 0 yields 

<K*>y>o)  « o 

in  the  entire  xy- plane  with  the  exception  of  the  circle  S and  the 
strip  Z (on  which  * suffers  a discontinuity). 

Since  on  the  strip  2 both  condition  (1.10)  and  the  condition 
derived  from  equation  (l.ll)  must  be  satisfied 


IL = ■ (§L.o 


and 


(tx)  " (fj)  “ 0 for  |y|  < a » x2  + y2  > a2 ; x < 0 (1.12) 


z=+0 


z=*-0 


Finally,  since  the  fluid  far  ahead  of  the  wing  is  assumed  to  be 
undisturbed,  the  condition  at  infinity  is 


lim 

X~N-» 


^5 

dx 


lim 

x-*+* 


i?  _ 

ay " 


lim 

x-h-« 


(1.13) 


In  the  hydrodynamic-  problem  under  consideration,  account  is  taken 
of  the  distribution  of  the  vortices  along  the  surface  of  the  wing.  It 
is  this  circumstance  which  makes  the  treatment  more  accurate  than 
the  usual  wing  theory. 


The  hydrodynamic  problem  is  thus  reduced  to  the  following  mathe- 
matical problem:  ‘ To  find  a harmonic  function  <p(x,y,z)  regular  over 
the  entire  half-space  z > 0,  which  on  the  circle  S satisfies  the 
condition 


(1.14) 
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on  the  strip  Z , the  condition 


on  the  entire  remaining  part  of  the  xy- plane,  the  condition 


(1.15) 


.(*>y>o)  - o (i.i6) 

and  the  derivative  a of  which  remain  hounded  in  the  neighborhood  of  the 
rear  aemicircumference  BCD,  while  in  the  neighborhood  of  the  forward 
semiclrcuoference  RAD  they  may  approach  infinity  as  where  3 

is  the  distance  of  e point  to  the  semicircumference  BAD.  Finally 
the  condition  at  infinity  (1.13)  must  be  satisfied. 

An  expression  for  the  harmonic  function  <p(x,y,z)  is  given  in 
closed  form  depending  on  an  arbitrary  function  f(x,y)  ~ satisfying  all 
the  imposed  requirements  besides  equation  (1.14).  The  function  £(x,y) 
can  be  determined  from  this  condition,  that  is,  the  shape  of  the  wing 
corresponding  to  the  function  f(x,y).  An  integral  equation  is  also 
given,  the  solution  of  which  is  reduced  to  the  determination  of  the 
function  f(x,y)  for  the  given  shape  of  the  wing,  that  is,  for  a 
given  function  C(x,y). 


2.  Derivation  of  the  Fundamental  Equation 

Inside  the  circle  ABCD,  the  point  Q with  coordinates  £,  q is 
taken  and  the  function  K(x,y,z,£,Tj)  constructed,  where  x,y,z  are 
the  coordinates  cf  the  point  P,  according  to  the  following  conditions: 

(1)  The  function  K,  considered  as  a function  of  the  point  P,  is 
a harmonic  function  outside  the  circle  ABCD. 

(2)  The  function  K becomes  zero  at  the  points  of  the  plane  xy 
lying  outside  the  circle  ABCD. 

(3)  The  derivative  dK/dz  becomes  zero  at  all  points  of  the  circle 
ABCD,  except  the  point  Q. 

(4)  When  the  point  P approaches  the  point  Q,  remaining  in  the 
upper  half-space  z > 0,  the  function  K increases  to  infinity  but 
the  difference  K - (l/r),  where 


r =*  O2  + (y  - n)2  + z2 


remains  bounded. 
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(5)  The  function  K remains  finite  and  continuous  In  the  neigh- 
borhood of  the  contour  C of  the  circle  ABCD. 

Because  of  the  second  condition,  the  values  of  the  function  1C  at 
two  points  situated  symmetrically  with  respect  to  the  plane  xy  differ 
only  in  sign: 


K(x,y,-z£,n)  - - K(x,y,z,C,n) 


(2.1) 


as  follows  from  the  principle  of  analytic  continuation.  It  is  then 
evident  that  if  the  third  condition  is  satisfied  on  the  upper  side  of 
the  circle  ABCD  it  will  be  satisfied  also  on  the  lower  side,  since 
according  to  equation  (2.l)  the  derivative  dK/dz  has  the  same  value 
at  two  points  situated  symmetrically  with  respect  to  the  xy- plane. 

It  is  evident  further  that  when  the  point  P approaches  the  ijoint  Q 
from  below  eo  that  z<  0 then  K(x,y,z,£,t|)  will  behave  as  - l^r. 

Because  of  the*  third  condition,  the  function  K can  be  continued 
into  the  lower  half- space  through  the  upper  side  of  the  circle  ABCD 
as  an  even  function  of  z.  Thus  a second  branch  of  the  function  K is 
assumed,  again  determined  over  all  the  space  outside  the  circle  ABCD 
and  differing  only  in  sign  from  the  initial  branch  of  the  function  K. 

It  is  then  evident,  however,  that  at  the  points  of  the  upper  side,  of 
the  circle  ABCD,  the  values  of  the  second  branch  of  the  function  K 
and  its  derivatives  coincide  with  the  values  of  the  first  branch  of 
the  function  K and  its"  derivatives  at  the  points  of  the  lower  side 
of  the  circle  ABCD.  That  is,  in  the  analytic  continuation  of  the 
second  branch  of  the  function  K through  the  upper  side  of  the  circle 
ABCD  into  the  lower  half-space,  the  initial  branch  of  this  function 
is  again  obtained. 

A two-sheet  Riemann  space  is  considered  for  which  the  branching 
line  is  the  circumference  ABCD.  In  this  space  K(x,y,z,£,T|)  is  a 
single-valued  harmonic  function  remaining  finite  everywhere  with  the 
exception  of  the  two  points  Q having  the  same  coordinates 
but  belonging  to  two  different  sheets  of  space j in  one-sheet  the 
function  K be  ’ /es  near  the  point  Q as  l/r  and  in  the  other  sheet 
as  - l/r.  Such  a function  K(x,y,z,C,Tj}  can  readily  be  constructed 
by  the  method  of  Sommerfeld  (reference  2).  In  this  way  for  the  case  of 
a two-sheet  Riemann  space  having  as  branch  line  the  z-axis,  a harmonic 
function  V(p,<p,z)  (p,<P,z  being  the  cylindrical  coordinates  of 

the  point)  is  determined  which  is  single- valued  and  continuous  in  the 
entire  two-sheet  space  with  the  exception*  of  the  points  Q and  Q’ 
having  the  cylindrical  coordinates  (pW ,zf)  and 
where  near  the  point  Q the  function  V behaves  as  l/r  and  near 
the  point  Qf  as  - l/r,  where 


r «^p2  + pf2  - 2pp ' cos ((p  - <p’ ) + (z  - z')2 
r'  »^p2  + pf2  - 2pp'  cos (<p  + <pf)  + (z  - z')2 
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ThlB  function  V has  the  fora: 
V 

vhere 


- i Aarc  tan /HI  - A arc 

n |r  -x  r*  1 a - x ' f 


a m — - — J(p  + p*)2  + (*  - *')*  > ^ - cos  f- g * ; T'  - 

8Vpp” 

Setting,  in  particular, 

-xj  r -/^P2  + o'2  ♦ 2pp»  cos  * + (z  - «•)* 


cos 


t ♦ t' 


yields 
V 

or  finally 


2 J la  7 x . I a - t ! 2 T 

■ — <ax:c  tan/ arc  tanj >*  — arc  tan — , ms  ■■ 

jtr  | o - t \a+Tj  *r  y o2  - * 


2 Z-V/PP1”  sin  f 

V ■ —arc  tan 

itr  r 


An  inversion  with  respect  to  the  point  with  coordinates  p - a, 
<p  - 0,  z « 0 is  carried  out. 

2a2(x,  - a)  Za2(Ci  - *) 

p cos  <P  ® a + - — ^ — ; - P = a + 


(xx  - a) 2 + yx2  + z ^ 


p sin  <P  » 


2a2yi 


(«1  - a)2  + Ci2 

2a^zn 


5 z « 


(xx  - a)2  + yx2  + z^2  (x-,  - a)2  + y^2  + z^ 


2a2  ^ 


(«!  - ®)Z  + 


The  function 
V,  = 


2aZV 


^ - a)2  + yx2  + If  ]/ ~ a)2  + £,* 


to 
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expressed  in  the  variables  x^,  is  then,  as  is  known,  a harmonic 
function*  Computing  it  and  replacing  _x1#  y^,  z±  by  y,  z,  x and  Ki, 
by  q,  £ yield  the  required  expression  of  the  function  K(x,y,z,£,q) : 


2 4/a2  - £2  - n2 

K(x,y,z,€,q)  » — arc  tan-* *- 

nr 


V a2  - x2  - y2  - 


z2  + R 


(2-2) 


ar 


valid  for  z > 0,  vhere 


r - V(*~-  O2  + (y  - *02  ” *2 

(2.3) 

R «^(a2  - x2  - y2  - z2)2  + 4a2 z2  «Y(a2  + x2  + y2  + z2)2  - 4a2 (x2  + y2) 

That  this  function  satisfies  all  the  above  set  requirements  is 
easily  verified*  the  arc  tangents  must  be  taken  between  0 and  n/2* 
for  z < 0 the  value  of  the  function  K is  obtained  by  equation  (2.1). 

The  following  function  is  set  up: 

*l(x,y,z)  - ^ JT  KO,y,z,«,t|)ffc,n)d«an  (2.4) 

s 


where  f(x,y)  is  an  arbitrary  function,  which  is  continuous  together 
with  its  partial  derivatives  of  the  first  and  second  order  in  the  entire 
circle  S,  and  the  integration  extends  over  the  entire  area  of  the 
circle  S.  Evidently,  $i(x,y,z)  is  a harmonic  function  in  the  entire 

space  outside  the  circle  S.  Because  of  the  first  property  of  the 
function  K,  the  function  (p^x^z)  becomes  zero  at  all  points  of 
the  plane  xy  which  are  outside  the  circle  S.  Hence  equations  (1.15) 
and  (1.16),  which  must  be  satisfied  by  the  solution  $(x,y,z)  of  the 
problem  posed  in  section  1,  will  be  satisfied  for  the  function  $j(x,y,z). 

The  function  <p^(x,y,z)  does  not  in  general  satisfy  the  condition  of 
the  finiteness  of  the  derivatives  of  this  function  on  -the  rear  half  of 
the  contour*  of  the  circle  S.  For  this  reason,  a function  such  that 
the  obtaii  ed  function  <p(x,y,z)  also  satisfies  this  condition  is  added 
to  <P1(x,y,z). 

The  following  equation  is  evident: 


= _i_  p p 

dx-  ” 2rt  J j 

b 


fj-  ie-,n)d£dn 

dx 


} 


! 
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The  character  of  the  approach  of  the  function  dK/dx  to  infinity  is 
considered  as  a point  approaches  the  contour  C of  the  circle  S.  As 
may  be  easily  computed 


If  a point  with  coordinates  x,y,z  is  near  the  contour  C of  the 
circle  S the  distance  of  this  point  to  the  contour  C is  denoted 
by  5 1 then 

6 a2  + x2  + y2  + z2  - 2a  Vx2  + y2  (2.6) 

Hence  near  the  contour  C,  the  approximate  equation  holds: 

R - 2a5  (2.7) 

When  the  fixed  point  £,t|  lies  inside  the  circle  3 while  the 
point  with  coordinates  x,y,z  lies  near  the  contour  C of  the  circle, 
then,  as  follows  from  equation  (2.5), 

^ a2  - C2  - T)2  Va2  - x2  - y2  - z2  + R + 0(l)  (2.8) 

«x  «ar^R 

where  the  symbol  0(l)  denotes  a magnitude  which  remains  finite  when 
6 approaches  zero.  Thus  dK/dx  has  the  order  l/V&*  The  principal 
part  of  dK/dx  is  not  a harmonic  function.  It  is  not  difficult,  however, 
to  find  a harmonic  function  having  the  same  infinite  part  near  the 
contour  C as  dK /Sx.  For  this,  it  is  sufficient  to  forn,  after  the 
analogy  of  equation  (2.5),  the  derivative  dK/d5j  this  derivative 
remains  finite  near  the  contour  C of  the  wing*  moreover  it  is  easy 
to  see  that 


\ 
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This  function  is  harmonic  and  differs  from  dK/dx  by  a quantity 
which  remains  finite  near  the  contour  C. 

By  computation,  it  is  further  shown  that  the  function  Just  described 
is  represented  in  the  form  of  the  ?:.tegral 


a*  ^ a*  ^ l P*  V*2  - c2  -"n2  V«2 

^ * V2  -*  (x2 4-  y2 ♦ *2 ♦ *2 - 2ax  oo«  r-2»y  «la  r)  (C*  >**♦**•  t*<  eo«  r-2»n  »m  r) 


where  the  function 


Va2  - x2  - y2  - 22  4-  R 

x2  + y2  + z2  + a2  - 2ax  cos  r - 2ay  sin  T 

is  a solution  of  the  equation  of  Laplace  having  the  circumference  C as 
the  branching  line  and  the  point  with  coordinates  (a  cos  y,  & sin  y,  0) 
as  a singular’  point.  From  this  it  follows  that  the  function 


5*  

V*2  - *2  - n2  V>2  • i2  ■ y2  » i2  ♦ r cot  r ir 
(x2  + y2  <•  e2  ♦ a2  - 2mx  cos  r - 2 ay  sin  r)(C2  ♦ n2  * »2  - 2*C  cos  r - 2»ij  sin  y) 

ajc  l *A2  ■ C2  - ,n2  V*2  - X2  - y2  - »2  + r co«  r dr 

a*  *2^S  J (x2  ♦ y2  ♦ x2  ♦ a2  - 2ax  cor  r - 2ay  sin  y)(*2  ♦ n2  + a2  - 2a*  cos  r - "in  r) 


(2.11) 


remains  finite  near  tbe  points  of  the  rear  semi circumference  of  the 
* circle  S. 

Therefore  it  is  assumed 


■ a» 


3k 

i Z7- iH  - 1 - n* 

k2V2 J (x2  + jr2  + t2  + »z  - cax  coa  r ■ 


_ . rr  L 

S * r.  JJ  '''’'ijs' 


Va2  - x2  - y2  - t2  ♦ R cob  r dr 
2ny  Bin  y)(*2  + H2  ♦ s2  - 2a*  cos  T ^2xij  tin  r) 


h 


dtj  (2.12) 
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Integrating  with  respect  to  x and  considering  the  condition  at  infinity 
(1.13)  yield  the  final  equation 


♦(«* 


[*(*»/!*#*#*)  ♦ 


% 3 *\ 

i f*  r v»8„-  ahm nnzi in » ~v  * * L**,^.*, 

J (*2  ♦ y2  ♦ f2  ♦ a*  - Z*x  cos  r - 2»y  sin  r)(&  ♦ n2  ♦ a*  - 2t(  co«  r - Sin  r)J 


This  equation  may  he  written  in  somewhat  different  form.  Because  of 
equation  (2.11)  — 


fi-  -i  «-=k»- 


ffj t 

8 '2 


x , ^LLz±^IZI^^l2Ll± r d, 

(kc  ♦ jr  + - £*x  COB  r - 2ay  tin  r)(C^  ♦ 1**  ♦ - 2ajr  co»  r - Ear)  tin  r) 


Since  the  function  K becomes  zero  on  the  contour  C 

JY  n *«.*>*<  *>-- ff  t2-«> 

s s 

Introduction  of  further  notations 

. 1 f f ^ZZEZZ g&atei *» G(>);  F(«,n) 

2iP^/2  v £2  + n2  + a2  - 2a £ cos  y - 2aTj  sin  y 95 


results  in 


(2.15) 


V&2  - x2  - y2  - z2  + R G(r)  cos  X 


r . 

J n x2  + y2  + z2  + a2  - 2 ax  cos  y - 2ay  sin  y 
*2 


dr  (2.16) 


g 
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and  after  integration  with  respect  to  x 


x 

K(x,y,z,g,n)dx  PU,»l)d^  dt) 

C 4-  • 


ft 


i 


.2  . ,2 


y2  - z2  -i-  R G(r)  cos  r 


+ y2  + z2  + a2  - 2 ax  cos  r - 2ay  sin  r 


dx  dy  (2.17) 


The  given  functions  in  the  circle  S and  the  function 
G(r)  In  the  interval  (-*/2,  */2)  completely  determine  f(£,q),  co 
that  the  equations  (2.13)  and  (2.17)  are  equivalent. 


The  equation  $(x,y,z)  obtained  satisfies  the  conditions  imposed 
in  section  1. 


This  function  is  evidently  a harmonic  function  in  the  entire  space 
exterior  to  the  circle  S and  satisfies  the  conditions  at  infinity, 
equation  (l.l3)..  From  equation  (2.12)  it  follows,  that  in  the  plane  xy 
for  x2  + y2  > a2  the  condition  is  satisfied: 

(**)  =0 

and  from  equation  (2.13)  it  follows  that  ~~ 


<p(x,y,0)  = 0 

in  that  part  of  the  plane  xy  which  lies  outside  the  circle  S and 
the  stiip  2. 


It  remains  to  prove  the  finiteuess  of  the  first  derivatives  of 
the  function  <p(x,y,z)  at  the  points  of  the  rear  semi  circumference  C 
and  to  determine  the  behavior  of  these  derivatives  on  approaching  the 
points  of  the  forward  semi circumference  C. 

In  considering  the  neighborhood  of  the  real’  side  of  the  circum- 
ference C,  equation  (2.16)  may  be  used.  The  latter  shows  that  5<p/dx 
remains  continuous  at  the  points  of  the  rear  half  of  the  circumference  C 
and  becomes  zero  at  these  points. 

The  behavior  of  the  derivatives  with  respect  to  y and  z of  the 
following  function  is  considered: 
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*(x,y,i)  - J'J,K(x,y,£,£,T))F(£,ij)d$at|  (2.16) 


near  the  contour  C. 


0-  JY  | F(5’')d< 

s 


(2.19) 


Similarly  to  equation  (2.9), 


dK  dK 


Za-Jz  - x2  - y2  - 


zc  + R 


Sy  + 3tl  n^Za2!^  + (a2  - £2  - ij2)(*2  - x2  - y2  - z2  * r) 


i2  . r2  . „2 


and  similarly  to  equation  (2.14), 


- r - n* 


(2.20) 


where  this  part  of  the  integral  remains  finite  everywhere  and  on  the 
contour  C becomes  zero. 

In  order  to  evaluate  the  remaining  part  of  the  Integral  equa- 
tion (2. IS),  the  following  two  integrals  are  considered: 


'!<-)'//  *177 


IS  - E2 


(a2  - £2  - T)2)(a2  - x2  - y2  - *2  + R) 


- //; 


a2  - £2  - r,2  ^a2!*2  + (a2  - £2  - ij2)(a2  - x2  - y2  - z2  + H)J 


Both,  on  account  of  the  symmetry,  depend  only  on  *sjJ77  and  Z) 
hence  without  restricting  the  generality.,  it  may  be  assumed  that  y ■ 0, 
x > 0.  The  distance  5 of  a point  with  coordinates  (x,0,z)  is 
introduced  to  the  contour  C: 


6£t! 
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Since 


6 ■ Y(a  - x)2  + z2 


R ^ lx2  + z2  - a2.! 


the  following  relation  will  hold: 


Jl(x,0,i)«  / / rVai"  *2  -V  d*dTi 

J J 2a2  [(x  - £)2  + r^2  + z2l 


pDlar  coordinates  are  introduced 


£ ■ p cos  6 ; T)  ■ p sin  $ 


:,0,z)  < r 

Jo  Jo 


I j pVa2  ~ p2  dp  d* 

0 Jo  2a2  [p2  - 2px  cos  & +’  x2  + z2) 


Since 


- 2ex  008  * + x2  + *2  V(d2  + x2  -t-  z2)2  - 4p2  x2 


J^x.o.zX  4 f gjjEg  dg 

a Jo  a/'p2  ♦ x2  " z2)2  - 4p2  x2 


For  x a 


j!Uo,z)<  4 / - -piaL-  p2  dP  — 

‘ J0  V(P2  + X2)2  - 4p2  X2 


* oVa2  - 

*2  . x2  - 


- P2  dp  ^ jl  / 

-P2  W„ 


VT2  - p2  = a 
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While  for  x 44  a#  use  i®  made  of  the  Inequality 


R > a2  - x2  - z2 


btcin  p p 

8 1 


to  obtain 

ji(* 


VT2  -2  -2 


K - vr  *g*n 

a2  (x  - £)2  + tj2  + t2_  + (a2  - t2  - q2)(a2  - x2  - z2) 


1 J c Va2  - p2  dd  dp 

" 2 j J a4  - 2a2  xp  cos  M p2(x2  + z2) 


(2.23) 


'0  v0 

^2 


*1  -ir 

/(ja4  + p2(x2  - z2 )]  2 - 4a4x2p2 


0 YL°~  T ^ * ;-l  " A vo 

The  following  inequality  results: 


r£l*2JLfl  dP  x * f _ .g^P  _ 
_ “4-^2 


Jl(x>y..a)  < J 


(2.24) 


The  second  integral  is  considered.  As  before, 

«4  f 


J2(x,0,z)  4 -£ 


P dp 


L Va2  - p2  V(p2  »■  x2  + z2)2  4p2x2 


^0 


For  x ^ a 


4 f „ /JfL  f- -I- 

a Jo  V^2  - P2)  jn  + *)2  ♦ *21  [Cp  - ♦■Ir]  a J0  V‘2  - P2V(*  - x)2  + ,2  *2» 

For  x<  a an  inequality  of  the  type  in  equation  (2.23)  is  used: 


JgCx.O/z) * Jt 


P dp 


0 ^(a2  • P2)[a4  + p2(x2  + z2)  + 2a2xp][a2  - xp)2  + p2z2j 

r? 

^ I p 

•Jo  A/(a2  - p2)[(a2  - xp)n  + *2P2] 


* 1* 
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If  z - x and  therefore  z->/2#  then 

J2(x,0,z)  < 


* / dp 


,2  2a2  z 


*2  <_*2_ 


p“ ft^V2 

but  if  0 < z < ft  - x , and  therefore  6 < (a  - x)-y/Z',  then 

jt(«.o,.)«4/  — ~B  ^ <\f  — t-f 

* J0  (.2-xp)Va2-p2  • J (a-xW.2-^2  »2(»-x)  ft22 


The  following  approximation  is  obtained: 


J2(x,y,z)  $ 


i^V2 


where 


6 -^(ft  -V^TIJ)2  + z2 
Near  the  contour  C 


(2.25) 


(2.26) 


R » 2a&  (2.27) 

If  this  relation,  the  evident  inequality 

(a2  - x2  - y2  - z2|  ^ R 

and  the  obtained  approximations  are  used,  the  following  approximation 
is  obtained  from  equation  (2.20): 


It  is  evident  from  equations  (2.19)  and  (2.2l)  that  near  the 
contour  C 


dy 


(2.28) 


68 


NACA  TM  1324 


The  following  derivative  is  formed: 

I-  Mg  KM** 


But 

dK  _ 2z 
da  irr5 

where 


arc  tan  A + 


2 A r ^ + z(a2  + x2  + y2  + z2  - R)  1 
* 1 + A2  [ ^ rR(a2  - x2  - y2  - z2  + R)J 


w ^a2  » g2  - T|2Va2  - x2  - y2  - z2  + R 


ar-^2 

Hence  if 

|f(C,tj)  | < M 

then,  on  account  of  the  inequality 


A <1 


1 + A2  2 


for  z > 0 the  approximation  results: 


!$  -€  2M 
oz 


dr\  + 


2,V?  azM(a2  + x2  +y2  + z2  - R)  I Va2-  g2  - q2  dgdq 

sR  a/el2  - x2  - y2  - z2  + R 2a2r2  + 


(a2  - g2  - q2)(a2  - x2  - y2-  z2+  R) 


Noting  that 


ff 


A?  d£  dq  < 2* 


and  making  use  of  approximation  (2.24)  yield 
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N < 4«m  ♦ 2 yg  m ■*>*  + j + ** 

* • rV*2  - x2  - y2  - *2  + B 


Since  for  e > 0 


« * - (a2-  a2  - y2  - z2)  ^ Vr-  a^+xS+yS  + ag 

2m 


Va2-  x2-  y2-  *2  + R Vh2  - (a2  - x2  - y2  - z2)2 

hence 


< 4*M  + M (a2  + x2  + y2  + z2  - r)  - m2  + x2  + y2  + z2 
Ids  I aR 

Row  when  the  point  P(x,y,z)  is  near  the  contour  C,  then  be  cause  of 


R - 2aB  ; | x2  + y2  + z2  - a2  | < R 


there  is  obtained 


"° 


(2-29) 

Equation  (2.36)  ib  again  considered.  Since  the  derivatives 


— Ya2  - x2  - y2  - 72  + R 
dy 


. yVa2  - x2  - y2  - s2  » R. 


±Va2  . x2  . y2  . z2  7 a , ilgj±  ♦ y2  .t  gj^22 

dz  P\£T.  x2  - y2  - z2  + B 


)(2.30) 


* ~r  (a2  + x2  + y2  + z2  - R)  - a2  + x2  + y2  »•  z2  J 

ZaR  t ✓ 

have  near  the  contour  C the  order  it  is  clear  from  equa- 

tion (2.16)  and  the  obtained  equations  (2.28)  and  (2.29)  that  at  the 
points  of  the  rear  semi  circumference  of  C there  is  the  estimate 


(2.31) 
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But  is  Is  then  evident  that  the  derivatives  bp  /by  and  bp/bz  are 
finite  at  the  points  of  the  rear  eeml circumference  C. 

The  behavior  of  the  derivatives  of  the  function  * near  the  for- 
ward 8emlclrcumference  C can  readily  be  determined,  starting  from 
equations  (2.12)  and  (2.15). 

The  first  of  these  equations  may  be  written  in  the  form: 

Va2-  x2-  y2-  z2  + R Q(r)  cos  { 

— " ~~  """  1 ™ " ~ u™lfc  """ 

x2  + y2  + z2  +a2  - 2ax  cos  y - Bay  *in  r 

(2.32) 

But  on  the  one  hand,  the  estimate 


S 

holds  for  the  neighborhood  of  the  entire' contour  Cj  on* the  other  hand, 
on  the  forward  semicircumference  C,  the  second  integral  of  equation  (2.32) 
evidently  remains  finite.  Hence  for  the  forward  semi  circumference  C the 
first  of  the  estimates  is  obtained 


while  the  latter  two  of  these  estimates  are  obtained  in  a similar  manner 
from  equation  (2.13). 

In  this  manner  all  the  conditions  which  must  be  satisfied  by  the 
function  $(x,y,z)  are  satisfied. 

The  shape  of  wing  to  which  the  obtained  solution  corresponds  is 
explained.  By  equation  (1.14) 


9 


(2.34) 


Hence  it  is  necessary  to  find  the  value  b$/bz  in  the  plane  of  the 
circle  S.  Both  sides  of  equation  (2.13)  are  differentiated  with  respect 
to  z and  then  z set  = 0.  On  account  of  the  very  definition  of  the 
function  K, 
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l 


11m 

E-HO 


f(c>,i)dcd,i  ■ llB 


*(K>*)*Z  dTJ 


2 ft f(x,y) 

(2.35) 


Moreover,  on  account  of  equation  (2.30), 


Llm 

mo 


'Si 


^a2  - x2  - y2  - z2  + 


for  x2  + y2  < a2 
for  x2  ♦ y2  > a2 


If  this  Is  taken  Into  account, 


vherc 


(|)  * - f(**y)  + «(y) 


(2.36) 


vet 


VV  _ - t2  - n~  co>  at  »i 

? . ',»ji  oo«  r - z*r  *jb  r)(C*  ♦ n*  ♦ •*  - *»\  oo«  r 


*M  sin  t) 
(*•»*) 


For  the  function  t(x,y)  the  following  expression  is  found: 

x 

?(*,y)  •\S‘  f(x'y>dx  - ^ * + «i(y)  («•») 


where  6i(y)  is  an  arbitrary  function  of  y. 

Thus,  for  the  assumed  degree  of  approximation,  the  bending  of  the 
wing  in  the  transverse  direction  produces  no  effect  on 'the  form  of  the 
flow. 


It  is  assvjned  that  the  shape  of  the  wing  is  given,  that  is,  tbv 
function  £(x,y)  and  therefore  the  following  function-are  given: 

& = M(x,y) 


c 


(2.39) 
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From  equations  (2.34)  and  (2.36)  it  is  dear  that 

f(x,y)  - M(x,y)  + g(y)  (2.40) 

Substituting  this  value  in  equation  (2.37)  and  introducing  the 
notations 


. • C2  - n2'  w(c.n)  co»  r dr  d»  aids  

V*2  ♦ y*  - a2(**>  y*  ♦ n1  - cob  r - 2*y  sla  r)(C*  ♦ **  ♦ •*"-  eoe  r - U*  ala  r) 

(1.41) 

* 

cot  y dr  d»  d| 

Yi*  ♦ y*  - »*(«*  ♦ y2  4 »*  • tu  cot  r • tin  y)(c2  ♦%*♦»*■  eon  t • M t) 


give  an  integral  Fredholm  equation  of  the  second  kind  for  the  determi- 
nation of  the  function  g(y): 

a 

g(y)  - N(y)  +J  Ii(y,n)g(n)dil  (2.42) 

-Q 


In  consideration  of  examples,  a function  f(x,y)  shall  be  given 
and  the  shape  of  the  wing  then  determined  by  equation  (2.38).  For  the 
obtained  shapes  of  the  wing  it  is  not  difficult  to  find  a solution  by 
the  usual  theory,  a fact  which  provides  the  possibility  of  evaluating 
the  degree  of  accuracy  of  the  usual  theory. 


3.  Computation  of  the  Forces  Acting  on  the  Wing 

The  fundamental  equation  determining  the  motion  of  the  type  under 
consideration  is  recalled: 


IT 

2 


V»2  - - H2V>2  • »2  - y2  - 1?  » R qob  r dr  dx 


(x2  4 yZ  + z2  + h?.  _ 2»x  cob  r - Zay  Bln  r)(t2  ♦ V'  «•  »2  - 2b£  cog  y - 2ai)  Bln  )r) 


(3.1) 
K(c,n)dCdn 
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The  vrJ.ue  of  the  function  9 for  the  joints  of  the  half-atrip  Z 
is  compute!.  Since  at  the  points  of  the  half-strip  Z 


this  value  is  a function  only  of  y.  The  notation  is  Introduced 

*(y)  ■ lim  *(x,y,z)  for  |yl  <a,  x2  + y2  > a2  f x<  0 (3.2) 
x-mO 

Then  evidently 

lim  f(x,y,z)  m - *(y)  for  |y|  < a,  x2  4-  y2  > a2,  x < 0 (3.3) 
r+-0 

The  circulation  over  the  contour  M'KM  (fig.  2)  connecting  the 
tvo  points  M and  M*  of  which  point  M*  lies  on  the  lower  and  point 

M the  upper  side  of  the  half-strip  £ , both  points  M and  M1  having 

the  sape  ccordinatee  x,y,0,  is  deno  ?d  by  T(y).  It  is  then  evident  that 


r(y)  »*(M)  -+(M')  - 2*(y) 


(3.4) 


Since  in  the  plane  xy  outside  the  circle  S both  the  function  K 
and  the  function  _ 


i2  - x2  - y2  - z2  + R 


become  zero,  it  is  clear  that 


m _i_  r P f*  P Va2  - c2  ~ - x?  - y2  r(c,n)  cob  r dr  dx  ag  dq 

2*3  J^2~2  Jn  (*2  y2  f a2  - Zax  cos  y - 2ay  flin  r)(£2  + H2  + *2  - T • 2*1 


2*q  sin  r) 
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Computation  shows  that 


f 5 


v»*  - E - E a ,/lH1*  *1bxL  .A  (s.*; 

♦ y*  ♦ a*  - 2u  cos  r - 2sy  sin  r If  1“  *ln  r ^ •' 


where  the  plus  sign  le  taken  for  y < a sin  r and  the  minus  sign  for 
y>  a sin  y. 

The  following  expression  Is  written  for  the  distribution  of  the 
circulation  In  the  vortex  layer  formed  behind  the  wing: 

3*  ~ 

*>■-$  [[?  - -A«.r *4*,  (3.7) 

* J J Jn  (&  ♦ I2  ♦ *2  - cos  r - 2*1  »m  r)  UM*  ,ln  r * y I J 

8 I 

The  forces  acting  on  the  wing  are  computed.  Denoting  by  p+  the 
pressure  at  a point  of  the  wing  S on  the  upper  side  of  the  wing  and 
by  pm  the  pressure  at  the  same  point  on  the  lower  side  gives  on  the 
basis"of  equation  (1.8) 


P-  - P+  * - 2pc 


(3.8) 


where  the  value  of  d<p/dx  is  taken  on  the  upper  side  of  the  wing. 
For  the  lift  force  P,  the  following  expression  is  obtained: 


P - j ' J1  ( P_  - p+)  dx  dy  « - 2pc  |?  dx  dy  - - 2pc  P dx  dy 

8 3 ^-a 

- - 2pc  [*(  V*2  - y^y^O)  - f(- V®2  - y2»y»oj]dy  - 2pc  J * f{y)  dy 


The  following  formula  is  obtained: 


P ■ pc 


a 


(y)dy 


(3.9) 


i 

V 


2439 


{ 

I 

i 


( 


I 


! 


i 


I 

I 

I 

i 


! 


NACA  TM  1324 


75 


having  the  same  form  as  in  the  usual  theory  of  a wing  of  finite  spe^. 

Bat  the  distribution  of  the  circulation  r(y)  by  the  present  theory  is 
somewhat  different  from  that  obtained  by  the  usual  theory*  The  deriva- 
tion given  is  not  connected  with  the  shape  of  the  wing. 


With  the  aid  of  equation  (3.6)  P may  be  directly  expressed  through 


3x 


. p . . SES£  f f f TTjrTtf  r(on)  f — - SSJLi'L, dt  d>i  (5.10) 

n2  J J K2  ♦ n2  ♦ a2  - 2a C eon  r - r 

8 f 

The  expression  for  the  Induced  resistance  W in  terms  of  the 
circulation  r(y)  likewise  has  the  same  form  as  In  the  usual  theory: 

w-feJ'  J*  r(y)  dy-dy'  <3-u> 

-a  -a  J J J 


because  the  origin  of  the  induced  resistance  is  due  to  the  fact  that 
behind  the  wing  a region  of  disturbed  motion  of  the  fluid  is  formedj 
the  kinetic  eneigy  of  this  disturbance  is  determined  on  the  other  hand 
exclusively  by  the  distribution  of  the  circulation  at; distant  points 
from  the  wing. 

The  expression  for  the  induced  resistance  is  obtained  from  the 
momentum  law. 

A surface  enclosing  the  wing  S is  denoted  by  Bj  the  momentum 
law  applied  to  the  wing  in  a steady  flow  then  leads  ter  the  expression 


«?  «■  f f pVnVxd<r  (3.12) 

B B 

where  n is  the  direction  of  the  outer  normal  to  the  surface  B and 
Vx,  Vy,  Vz  are  the  components  of  the  velocity  in  the  relative  motion 

of  the  fluid  about  the  wing.  Thus 


w”  ff  p cos(n,x)d 


- c + 


$2 

dx 


Vn  - - C 


cos(n,x)  + I? 


I 


t 

! 

i 


! 
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Substituting  these  values  in  the  preceding  formula  and  noting  that 

J ' J"1 cos(n,x)da  ■ 0 y J*  fa  m ^ 

BU  B 

results  in 

w " ' Iff  [(If)2  + ($$*  + (§)  + pff££da  (3,13) 

B B 

The  surface  B consists  of  a hemisphere  of  large  radius  with 
center  at  the  point  x - Xq  < -a  of  the  x-axis  enclosing  the  wing,  and 
of  the  circle  cut  out  by  this  hemisphere  on  the  plane  x « Xq.  With 
increase  in  the  radius  of  the  hemisphere  to  infinity  the  corresponding 
ports  of  the  integrals  entering  the  preceding  formula  approach  zero.  On 
the  surface  x * Xq 


' cos(n,x)  « - 1 ; 

on  ox 

therefore 

where  the  integration  extends  over  the  entire  plane  x * Xq.  For 
Xq  ■ «•  the  following  equation  is  obtained  : 

"-f/T[(l)2-(S)2]a'a*  (5-i5> 

where  $(y,z)  denotes  the  velocity  potential  of  the  plane -parallel 
flow  which  is  established  in  the  transverse  planes  far  behind  the  wing. 


The  usual  transformations  by  Green's  formula  yield 

a 

W =*  - P J*  $(y)  dy 
-a 


(3.16) 


where  the  integral  is  taken  over 
in  the  plane  yz. 

Since 

r (y)  = 2*(y) ; 

equation  (3.11)  is  obtained. 


the  upper  side  of  the  segment 


(-a,a) 


(3.17) 
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In  order  to  find  the  center  of  pressure,  the  principal  moments  of 
the  pressure  forces  about  the  Ox  and  Oy  axes  are  determined. 

For  the  moment  about  the  Ox  axis, 

a 

M*  - P f‘  (p.  - p+)y  dx.dy  - - 2pc  rpi  y dx  dy  - 2pe  J * (y)y  dy 

« " e -a 


from  which 


Mx  - pc  J*  yr{y)dy 


(3.18’; 


-a 


Expressing  Mx  in -terms  of  f(x,y)  yields 
5x 


M - - i ££ 2!  P P r2  - *2  - *i2  f(S,n)3in  r cob  f dr  dg  dn  (3  19j 

3 J J Jn  C2  + tj2  + a2  -_2a£  cos  r - 2arj  sin  y* 
s 2 

For  the  moment  about  the  Oy  axis, 

My  - - J1  J*  (p_  - p+)x  dx  dy  o 2pc  P J*  * dx  dy  (3.20) 

S 3 

Suocoituting  the  value  bf/bx  and  integrating  yield 

5k 

} V»2  - c2  - f(c»n)dc  dn  (3*21) 


cos2  r dp 

$2  + tj2  + . gaC  cor  r - 

2a»i  sin  r 
> 

The  following  values  are  obtained  for  the  coordinates  of  the  center 
of  pressure: 


- - ^ i yc  - ~ (3.22) 


4.  Examples 

NACA  comment:  Errors  in  the&e  examples  ere  lef erred  to  and  cor- 
rected in  the  paper  "Steady  Vibrations  of  Wing  of  Circular  Plan  Form". 
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The  equations  Just  obtained  are  presented  again: 

The  velocity  potential  for  z > 0 Is  determined  by  the  equation 


- ±JJ  ^ » 

- cob  y dr  dai 

c cob  x “ 2»y  Bln  r)(C2  ♦ t)2  + a2  - 2*4  00a  r - 2ai)  aln  r) 


(«.i) 

f(t,n)dt  «•) 


Kfx,y,z,5,T))  - —arc  tan  ^ ~ f *2  - £ ~ *2  + * 

" V2  « (4.2) 


R - Y(a2  - x2  - y2  - z2)2  + 4a2z2  * r =»  V (x  “ O2  + (y  - O2  + z2 

For  the  circulation  distribution  in  the  vortex  strip  formed  behind 
the  wing, 

rW  - - ^ » 

3* 

f r f V»2  n2  V*2..^.?..-.  y2  

J J JjfiCpJ*  (x2  + y2  + »2  - 2**  cod  r - 2ay  sin  r)(c'  ♦ I2  ♦ »2  - 2H  coa  r - Z*r\  aln  7) 


3* 

i r r r ^ 

* J J J*  (<?  + n2  ♦ *2  - cot  r - 28^] 


/all  1 alp 


aln  r)  I Vl  a *ln  r ‘ y I 


- 1 dr  dc  di| 


where  the  plus  sign  is  taken  for  y < a sin  r and  the  minus  sign  for 
y > a sin  r 

The  following  expression  gives  the  lift  force: 

3n  - 

P=pc  par(y)dy=-^  PPf2  ii.  - - 7 12  f(E>^8  T dr^dt|„ 

**  n «<,«  u , E’  + ij2  + - 2a£  cos  y - 2ai)  sin  y 


(4.4) 


The  usual  expression  for  the  induced  resistance  is 


a a 


w = ■§-  P P r (y)  ^ dy  dy’  (4.5) 

4n  J J dy’  y - y 


6£f<J 
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The  coordinates  of  the  center  of  pressure  are  determined  by  the 
equation?  - - 

*c  - - yc  - — - (4.6) 

P P 

where 


M,  - pc 


r*rT(r)V  . - 4 M»!  f f r **  • **  • **  r(c’l)“°  r”,r<r  llH  (4.7) 

J-a  3 *2  J J J«  ♦ V*  ♦ »2  - *•€  «»  r - 2*n  »io  r 


8 * 


.i£S 

* 


V. 


co«z  r dr 


a2  - z«(  coi  r >ib  r 


> V*2  - <2  - *)2  iln  (4.8) 


If  y is  set  equal  to  -a  cos  0 and  T (y)  is  represented  in  the 
form  of  a trigonometric  series. 


r(y)  “ Ai  Sin  0 + A2  sin  20  + . . . (0  < 0 < it)  (4.9) 

P,  W and  are  directly  expressed  in  terms  of  the  coefficients  of 
this  series  by  the  formulas 

— 

P - A1;  V = i„p£  rAn2  ; M*  = - i itpca2^  (4.10) 

8 n=l 

Finally,  the  shape  of  the  vine  is  determined  by  the  equation 

5(x>y)  " i f f(x,y)>ix  - x + 6l(y)  (4.n) 

c Jo 


where 


g(r) 


Vx2  + y2  „ a2  (j;2  + y2 


Va2  - g2  - n2  r(g,q)coc  r dr  Ox  dg  dq 
+ a2  - 2ax  ros  f - 2ay  sin  r)(g2  + n2  ♦ *2  - 


(4.12) 


Zag  coo  r - 2a rj  sin  j) 
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The  examples  are  now  considered. 

1.  First 

f(x,y)  - ca 


vhere  a Is  a small  constant. 

Polar  coordinates  are  used  and  the  following  integral  computed: 

f f- Va2  - e2  - r?  atdn 

u \ J £2  + 1]2  + a2  - 2a£  cos  r - 2ar)  ain  r 

m I Va2"-  P2  p dfl  dp  m f 2np  do 

Jo  Jo  a2  + pZ  • 2ap  TOs(*  * r)  JQ 

Substituting  this  value  in  equation  (i.3)  yields 

3* 


2it&  (4.13) 


r (y)  « - 2B£St  I cos  r 

ft 


VA  mm> 

f COS  r ./a(l  t 8lnr£  . x 

J«  |_V  | a sin  r - y | 


dr 


If  the  integral  is  taken, 

r(y) 


ca  J 

" V 


4a  + 2 ^a(a  - y)  f 2 ^/2~/a  + yj  - 

(a  + y)  log  ' ^rr~  - (a  - y)  log  \ (4.14) 

V2a  + va  - y V2a  t-A/a  + y J 

Setting  y = -a  cos  0 and  expanding  r(-o  cos  0)  In  a trigono- 
metric sine  series  in  the  interval  0 < 0 < * give  after  simple 
computations 


f(-a  cos  0) 


ai^^  4 + 


0 0 
4 cos  - f 4 sin  - 


(l  - cos  0)  log 


JL  ~ cos 


- 0 + COS  0)  log 


1 » sin 


1 + cos  ^ 
2 


1 + sin  - 


A1  sin  0 + ^ sin  36  + ^ sin  50  + . , . (0  * 9 < «)  (4.15) 


&> 


&v- 

MC 

v 1 
W« 

« v . 

‘tv ' 

%> 

ll* 

$ 

ijfe 

w- 
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where 


A 16a ca 

A1  * 2"  j A2k+1 


4aca  /l  1 1 \ 

[ - + - + ...  ♦ — = — j 

+ l)(2k  4 1)  \3  5 + 1' 


n2k(k  + l)(2k  4 1) 
(k  ■ 1#  2,  . . ) 


so  that 


16a  ca 


496a ca 
4725k2 ’ 


(4.16) 


I 

ill 

$i 

m 

& 

I 


1 

y 


The  distribution  of  the  circulation  obtained  is  very  near  that  of 
an  elliptical  distribution. 

The  lift  force  and  the  induced  drag  are  obtained  by  application  of 
equations  (4.10). 


P « ~ pcaAi  ■ — pa^c-a  * 2.5465  pa2c2a 

2 it 

W ■ ^ np(A^2  + 3Aj2  + . ..)  • 1.034  pa2c2a? 


(4.17) 


In  order  to  determine  the  position  of  the  center  of  pressure,  My 
must  be  computed  by  equation  (4.8). 


Equation  (4.13)  gives 


4 p 

My  = " 3 pcra°a  f 


My  Jt 

cc  " - P - £ 


(4.18) 


The  distance  from  the  center  of  pressure,  which  evidently  lies 
on  the  Ox  axis,  to  the  leading  edge  of  the  ving  thus  constitutes  about 
0.238  of  the  diameter  of  the  wing. 

In  order  to  determine  the  shape  of  the  wing  corresponding  to  the 
assumed  function,  it  is  necessary  to  form  the  function  g(y)  by  equa- 
tion (4.12).  If  equation  (4.13)  is  considered, 


/a2-y2  n 2 


, cos  y dy  dx  

fx 2 + y2  - a? (x2  + y2  + a2  - 2ax  cos  y • 2ay  sin  y) 

(4.19) 
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•llis  computation  shove  that  for  x > Va2  - y2 
3x 
pT 

I cos  r dr  _ . 


r — 

Jir  x2  + y2  + 


a2  - 2 ax  cos  r - 2ay  sin  r 2a(x2  + y2) 


log  ?l±.(y  -_»)2  + 


2a(x2  + y2)  x2  + (y  + a)2 


x(a2  + x2  + y2)  x2  + y2  - a2 

_i - arc  tan - 

l(x2  + y2)(x2  + y2  - a2)  2ax 


(4.20) 


y ■ - a cos 


9 l Va2  - y2  - a sin 


*o(e)  - — 8ln  M-a  cos  0) 


(4.21) 


for  0 < 0 < * 
✓ssi'n  0 


am  e J ttt  - cos  e t +4  c°«*  j 

t2-sin20 V 2(t2+cos20)~  2(t2+cos20)  1<>8  t2+4  sin4  | + 


t(t2  + 1 + cos20) 

(t2  + cos20)(t2  - sin20) 


■arc  tan 


t2  - sin20 


Computation  of  this  integral  gives 


Hq(0)  - ( - 2 


arc  tan  y 

Vi  - y2 


dy  ) sin  0 + 


l ( 1 + sin  / 1 + cos 

- Bin  01  log £ ] + ~ sin  0 [ log  — — & 1 + 

8 V !-■!»!/  8 V * -■»■{/ 

, . 0 .9 

1 - sin  — 1 - cos  — 

cos  | log  1 + sin  | log  1 

2 1 + sin  - 2 1 + cos  - 

2 2 
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The  shape  of  the  wing  Is  thus  determined  by  the  equation 

C(x,y)  - « [i  - ♦ ijf1  ^==  - 

_L  /i0g  VUi  2/Htf  - J_/0g  VE 1 ^EZj)2  - 
e*2  ^ Vs  - V^+y  8"2  \ VS  - -jrryJ 


2a  - ya  + y 


:y— log  TLT.,ufg  : v iog  ~v^a  " Va  " / 

2jt2V^  + y -^a  + Va  + y 2jr^-^a  - y -^5a  + -Ja  - y 


(4.23) 


This  wing  differs  little  from  a plane  wing  inclined  to  the  xy- plane 
by  a small  angle  a and  may  be  obtained  from  such  a plane  wing  by 
twisting.  The  values  of  the  function  £(x,y)  for  the  mean  value 
y ■ 0 and  for  the  values  y - ± a/2  are 


S(x,0)-«||  + £ / ^^dy.^log2(^+l)  + 


iVl 


log(V2  + 1)  - 0.8452  ax 


t («,  * I)  - ax  1 1 . 4 / dy  - i i„g*<2  , 'VS,  . -4  , 

' 2/  I 2 „*  / _ y2  2«2  - 8«2 

2 0 1 1 

-r log (2  + -^3)  + — log  3 j • 0.8335  ax 


It  is  of  interest  to  consider  what  results  for  the  obtained  wing 
are  given  by  the  usual  theory.  The  circulation  obtained  by  this  theory 
is  denoted  by  ^(y)*  if  the  expansion  of  this  circulation  in  a trigono- 
metric series  is 

r0(-a  cos  0)  = B-l  sin  Q + B2  sin  20  + , . . (5  < 9 < n)  (4.24) 

then  the  usual  theory  gives  an  equation  for  determining  the  coefficients 
which  in  the  case  considered  reduces -to  the  form 
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sin  n0  * 2icca  sin  0 


i ^ 


j _ sin  n0l 

U .BB"  "nr^4-25) 


Equation  (4.21)  yields 


t ^ * f) 


. . 

sin  n0  ■ 2ncaa  sin  9 • — — Bq(0) 


(4.26) 


Expansion  of  the  function  Hq(0)  into  a trigonometric  series  is 
sufficient  to  determine  the  coefficients  Bq.  Despite  the  complicated 
form  of  the  function  Hq(0),  it  can  be  expanded  and  in  the  interval 
0*  0*  n / * - v 

* _ . I /i  rtp  tan  v 1 


Ho(e)  - sin  0 M . 4 - 2 f dy  + 

\16  J0  VT77  ) 

> ein(2k  + 1)9  T , . 1 . 1 2(2k  + l)2  + 1 1. 

k(k  + 1)  ['  3 *•*  + ITT I ‘ (4k  7 l)(4k  + 3)J^4*27) 


Ho(e)  - Lt  P2!t+i  sln(2k  + i)® 

to=0 


» - 0.1389  \ P3  = - G \ P5  » - 0.1213 
07  * - 0.0460  j 09  = - 0.0212,  ... 

Equation  (4.26)  shows  that 

4aca(*2  - Pi)  4acaP2k+l  / « % 

Bl  ■ J ®2k  a 0 I = “ ” (k  * 2>  • • • ) 

*(*  + 2)  2 + n(2k  + 1) 

(4.20) 

The  numerical  values  of  the  first  coefficients  will  be 


= 2.4784  oca  ; Bj  = 0.0562  oca  > B5  « 0.0087  oca 
B7  ■ 0.0024  oca  I Bq  *=  0.0009  oca,  ... 


6ff»2 
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The  following  value  is  obtained  for  the  lift  force  2 


Pq  « ^ JtpcaB^  « 3.8932  pc2a2a 


exceeding  the  accurate  value  by  53  percent. 


FOr  the  induced  drag. 


Wn  * 2.416  pa2c2a2 


with  an  error  of  134  percent. 


2.  If  a is  assumed  to  be  small,  f(x,y)  * - 2cax  is  taken. 


The  circulation  r(y)  is  computed.  First  the  value  of  the  fol- 
lowing integral  is  found. 


u g u z*  + Vr  + a2  - 2a£  cos  r - 2arj  sin  r 3 


i2  - r2  - n2 


Equation  (4.3)  gives 


8ca2a 


cos2r  \l-^1  i6inr)  .-1  dr 
_V  |a  sin  r - y| 


The  computation  of  this  integral  leads  to  the  very  simple  expression 


r(y)  = 2ca(az  - y2) 


(4.32) 


Thus  in  the  case  considered,  a parabolic  distribution  ox"  the  circu- 
lation was  obtained.  For  this  reason  the  computation  of  the  forces  can 
be  easily  carried  out: 


na  e 

P *»  pc  I r(y)dy  = - apc2a3  • 2.667  apc2a3 


(4.33) 


W - - pc2a2a4  * 1.2732  oc2a4a? 
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Equation  (4.3l)  Is  used  In  the  computation  of  My  by  equation  (4,8): 

My  - — pci 2a4a«  1.509  pe2a4a  » xc--liK--iSa  (4.34) 

1 27*  c P 9k  ' 

In  order  to  de  to  mine  the  shape  of  the  wing  it  is  necessary  to 
compute  the  function  g(y)j  equation  (4.12)  yields 


fS- 


5 pVa^-r  pZ 

J+-  J« 


cos*r  dy  dx 

y2  - a2(x2 + y2  -fa2  - 2ax  cos  y-  2ay  sin  r) 


Setting 


Hi  (6)  2-  sin  8g(-a  cos  8)  (0  < 8 ^ *) 

x 4dca 

and  carrying  out  the  integration  with  respect  to  r yield 


(4.35) 


Hi(e)  - 


sin  0 


(t2  + cos28)2  Vt2  “ sin28 


(t(t. 


2 + cob28)  + 


i t2  + 4 0084  *9 

- (co828  - t2)(t2+  1+ cob28)  - - t COB  8(fc2+  .1+  C0826)l0g  — — - \ 4 

1 2 t2+  4 Sin4  £ 

2 


2(t2  -f  co828)  -f  (t2  - COB  2e)[i  + (t2  + cos  fafl 
2(t2  - tin2©) 


arc  tan 


t2  - sin28\ 

2t  y 


Integration  yields 

r 


«(i - !)  * i w * 

^ V J V ‘ c°81)  * 8ln  |) 


sin  8 cos  8 log  lan-  + - log 


(l  + cos  |)  (l  - 


(l  - cos  ?)  (l  + sin  |) 


(4.36) 


In  equation  (4.11),  the  following  is  taken: 


- gi(y)  = a(a2  - y2) 
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Then  for  the  function  £(x,y),  which  detemlnec  the  shape  of  the  wing, 
the  following  expression  is  obtained: 


C(x,jr)  - a(a‘ 


y2) 


2oax 

x 


11 

y) 


L \IEZ2  \Fry , -1-  l0,  1:  ^ V* 

\ V 2a  V 2a  12  8 (^w<1  . JT^y)(-jX  - 

x logjm . j.  io8  (Vs-*  .v^zy)(VF-j^m\  (4>s. 

a Va  - y 4a  (^/JJa  - Va  - y)(V^*  ♦ ~Va  + /)  J 
This  wing  is  thus  obtained  as  a deformation  of  the  wing: 

C(x,y)  - a(a2  - x2  - y2) 

which  for  small  a differs  little  from  a segment  of  a sphere. 

In  particular,  for  y » 0, 

5(x,0)  - a(a2-  x2)  + [1-V2+5  log(V2+l)]-  a(a2- x2  - 0.0767ax) 

In  order  to  apply  the  general  theory  to  the  obtained  wing 
is 'expanded  into  a trigonometric  series: 

Hi(0)  - ^ ^ - J 3 

£ f.  12nk(k  + 1)  + 

teel  4k(k  + l)(2k  - l)(2k  + 3)  |_ 

2(16k2  + 16k  - 3)  ^1  - | + i - ...  + -—jj*  6(2k  + 1)J 


log  tan  7;  1 ®in  ® + 


£ 

k*0 


rzk+isln(2k  + x)0 


(4.38) 


where 


r1  = - o.693i  1 r3  * - 0.1783  1 r5  - - 0.0812 
Xy  - “ 0 . 0463  j Yg  » - 0 . 0300 , ... 
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For  tin  case  considered,  the  usual  theory  gives  for  the  determination 
of  the  circulation 


r0(-a  cos  8)  m Bj  sin  8 + Bg  sin  28  + ...  (0  ^ 8 < x) 


the  equation 


Bjj  sin  n8  ■ 2icca  sin  8 \aa  sin  8 


8 (-*  COB  0)  1 > Bin  n*l 

; isfe"s»-nr5j 


(4.39) 


Equation  (4.35)  and 


sin20  .-51  6ln(2k  4 1)9 (o  « 6 < it)  (4.40) 

" feo  (2k  - l)(2k  + l)(2k  + 3) 

give  from  equation  (4.39)  the  equation 

m m 

S + t)  = - (2k  -~)(2Trf)T2>r"  3)]8ln^2k  + v0 


(4.41) 


from  which  without  difficulty  ^ is  obtained,  in  particular 
B2k  = 0 j Bi  = 1.8457  aca2  ; B3  = - 0.2132  oca2 
B5  = - 0.0250  oca2  j B7  = - 0.0075  oca2  ; Bg  «■  - 0.0032  oca2. 


The  following  value  is  obtained  for  the  lift  force: 


P =*  npeaBj  = 2.899  ac2a3p 
exceeding  the  accurate  value  by  8.7  percer.  . 


(4.42) 


The  induced  drag 


W = 1.3927  pa2c2a4 


(4.43) 


exceeds  the  accurate  value  by  9.4  perce  it. 
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where 

. 128  asa2 

. wf — sdii — (».i. ... 

*2  lsk(k2  - l) (4k2  - 1)  \ 5 


(4.47) 

1 ^ -2*  1 
4k  - 1/  (k2  - l){4k2  - 1)J 


* 2|3;  •••) 


80  that 


Ag  *»  * 0-4803  oca2  j A4  ■ 0.00549  aca2 


*»  0.00234  ace2  i Aq  - 0.00121  aca2 

Evidently  there  ic  hr  lift  force,  whereas  for  the_ induced  dre.g 
the  following  value  is  obtained: 

V - o .1813  pa2c2a4  ~ (4.48) 

The  moment  cf  the  forces  about  the  Ox  axle  is: 

r *■  t « 0.3V iZ  (4.49) 

‘fhe  moment  of  the  forces  atci<t  the  Oy  axis  io  computed  with  the 
aid  of  equation  (4.8)  0 where  use  miiyt  be  mp.de  cf  the  result  (4.d4), 
it  is  found  that 

My  * r (4.50) 

The  fa  llowing  function  it*  ’xov  cuurnu  .ed: 

co 3 y dy  dx 

+ a2  - 2ax  cos  y - 2ay  sin  y) 
(0 <9-?  v;)  (4.51) 


e(y; 


Pr^rjc 

3,2 


_ sin  r 
Vx2-*  y2-e2(x2  + y2 


Setting 


n2(e) 


rin  C^a  cog  9) 

2 s vu? 


V '.i.  ■••  » t.&  -fctfjfc.  j-  ■ _j?  ‘ .g-  :V*  -■•  . •:  >-  ■•  - 
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and  carrying  out  the  integration  vtth  respect  i*,  y give 


*nin  0 


a?.(o) 


sin  6 


(t2  + cpb20)2  Vt2  - sin2#) 


cob  tf(t2  + cob20)  + 


i t C^,3  0(t2  + 1 •*■  CO*I20)  + 


j (t2  + 1 + COB20)(cOS20  - t2)  log  — 


t2  + 4 COS*  | 


t co-  0 [l  + (t2  + COB  20}2] 


t2  - ain20 


—arc  tan 


2 + 4 Bin*  1 
2 


t2  ^ aln20\ 
Zt  J 


sin  0 con  0 


(4.52) 


(f  . O'? 

1 + COB  - \ 

log 2 j 

v »-«!/ 


1 + 3 cos  0 . 0 . 1 + 005  2 

— * — 2 108  rrr® 


(.  0>2 

log  1 1 Bln  I ) + 
> 1--!/ 


1 - 3 COB  0 0 1 + Bin  - 

cos  - log 

2 2 1 - Bin  * 


Expansion  ^ a tricono.ee trio  serlec  gives 


n2(0)  »■  sin  20  J log 2(V2  + i)  + 


3 / arc  tan 

u 


y 9*2  2159 

fdy"T(T+  630~ 


) 2(dk2  + 1)  /.  i 1 12k2  \ , 

Lj  — -i ~L ( 1 + t + . ••  + T2 T - — 3 j sin  2k0 

krj  (f^2  . ;)(4k’  . 4)  V 3 4k  - 1 8k2  + 1/ 


52!c  sin  2kfl 


(4.53) 


r TyyyTvww^  v rmrn  ;zwr  vrei  ^urrrrr/v*mrrv*vin‘*\’  -*  ^rrvTZ'TV  / dcwtgvji 


NACA  TM  132:4 


where 


&2  « - 0.27412  | S4  - - 0.06127  ; &6  - - 0.0519Q  j 5Q  « - 0.03641, 


The  usual  theory  for  determining  the  circulation 


m 

r0(-a  cos  0)  - iC  Bn 


sin  nd 


gives  the  equation 


ct 

Yj  BnO-  + t) 


tin  nd  ■ 2nca  sin  6 <-a a cos  0 - 


■a  cos  6)1 
c / 


o 4uca“ 

■ - ncaza  sin  26  - — HgC©) 


(4.54) 


from  which  without  difficulty 


Bofc+i  * 0 (k  = 0,1,2,  ...) 


B2»-  0.7304  aca2  j -0.0047  aca2  j Bg«  0.0021  uca2  r Bg  *0.0011  aca2, 


The  lift  force  is  found  equal  to  zero  and  the  induced  drag  and 
moiaent  of  the  forces  about  the  Ox  axis  are 


W - 0.4191  pa2c2a4  ; K,.  0.5737  pa2c2a4 


(4.55) 


The  first  givee  an  error  of  131  percent,  the  second  of  52  percent. 


3/  a combination  of  the  obtained  solutions  it  would  have  been 
possible  to  obtain  further  examples.  From  the  examples  given  it  Is 
clear  that  for  the  case  of  a circular  wing  considerable  devis.^ions  are 
obtained  vetween  the  usual  ard  tne  exact  theories - 
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